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Abstract 

We study a gaseous atomic Bose- Einstein condensate with laser-induced dipole-dipole interactions 
using the Hartee-Fock-Bogoliubov theory within the Popov approximation. The dipolar interactions 
introduce long-range atom-atom correlations which manifest themselves as increased depletion at 
momenta similar to that of the laser wavelength, as well as a 'roton' dip in the excitation spectrum. 
Suprisingly, the roton dip and the corresponding peak in the depletion are enhanced by raising the 
temperature above absolute zero. 
PACS: 03.75.Fi, 05.30.Jp 

1 Introduction 

One of the novel features of gaseous Bose-Einstein condensates (BECs), when considered from a many- 
body physics perspective, is the ability to directly manipulate the interatomic interaction using external 
electromagnetic fields pQ. In the BECs realized thus far the interatomic interactions can be described 
by a single parameter, the s-wave scattering length, a, and are short-range in comparison to the average 
interatomic distance. The experiment by Inouye et al demonstrated how the s-wave scattering length 
can be changed by magnetic fields via a Feshbach resonance. Here we consider enhancing and controlling 
the interatomic interactions using laser-induced dipole-dipole forces. These interactions are intrinsically 
long-range and so affect the gas in a way profoundly different from a Feshbach resonance. In particular, 
dipole-dipole interactions, both laser-induced PJ and static 0], have been predicted to be capable of 
introducing a 'roton' dip into the Bogoliubov excitation spectrum of the BEC. This is behaviour remi- 
niscent of that of the strongly correlated quantum liquid helium II, and can be understood in terms of 
Feynman's relation which provides an upper bound for the spectrum E(k) of helium II 

where S(k), known as the static structure factor, is the Fourier transform of the pair distribution function. 
S(k) is a measure of 2nd order correlation in the system. Feynman's formula interprets the roton dip in 
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helium II as being due to a peak in S(k) and hence as being due to strong correlations, which are at values 
of k corresponding to phonons having wavelengths on the order of the interatomic spacing (intuitively: the 
excitations at these wavelengths require less energy to excite in comparison to neighbouring wavelengths 
which cost more compressional energy). 

In a gaseous BEC at zero temperature the Bogoliubov theory of a weakly interacting degenerate Bose 
gas [5] is valid and shows that 2nd order correlation arises predominantly from pairs of atoms scattering 
out of the condensate to form the so-called quantum depletion, or 'above condensate' fraction. In a 
homogenous system the atoms in these pairs are prefectly correlated, in the sense of perfect two mode- 
squeezing |H1[7||H]. However, in gaseous BECs the quantum depletion, and hence 2nd order correlation, 
is typically very small on account of their diluteness, as expressed through the gas parameter no? -C 1, 
where n is the atomic density. Indeed, within the Bogoliubov theory the Feynman relation Q is an exact 
equality, yet in a measurement of the bulk excitation spectrum of a regular gaseous BEC interacting via 
s-wave scattering by Steinhauer et al jS] no roton was seen. Laser-induced dipole-dipole interactions, on 
the other hand, introduce controllable (via laser intensity etc.) long-range correlations on the scale of the 
laser wavelength, and one might expect some significant depletion at the corresponding momentum. 

In our previous investigations into laser-induced dipole-dipole interactions in a BEC (and 
references therein) we considered the zero-temperature case, and the calculations were limited to the 
basic Bogoliubov method [5] which does not treat the depletion self-consistently. Our purpose here is to 
use the so-called Popov version of the Hartree-Fock-Bogoliubov method [1 11112) , which treats the depletion 
as a self-consistent mean-field in order to include the effects of the the back-action of the depleted fraction 
upon the condensate. This method was primarily developed to treat the finite temperature case where 
depletion is important, but our hope here is to account for significant quantum depletion due to enhanced 
interactions and also for thermal depletion due to non-zero temperature. The most significant results 
we find here are that, perhaps contrary to expectation, the depletion can increase the effects of the 
dipole-dipole interactions, deepening the roton dip in the spectrum, so that the roton is not diminished 
by finite-temperature effects but can actually be enhanced. 

2 Laser induced dipole-dipole interactions 

The idea of exploiting dipole-dipole forces to modify the properties of atomic BECs first arose in the 
context of static dipolar interactions. A BEC of atoms with a large permanent magnetic moment (e.g. 
chromium) ordered by an external magnetic field, or equivalently, polarizable atoms in a static electric 
field, has been considered by a number of authors EI El EH El EE] • Here, however, we are interested 
in fully retarded dynamic dipole-dipole interactions, such as those induced by an electromagnetic wave 
(e.g. laser beam). The dynamic dipole-dipole interaction is distinguished from the static case by a longer 
range (the retardation gives it an attractive r~ 1 component which can be used in certain geometries to 
mimic gravity |19p and a huge enhancement of atomic polarizabilty close to a resonance. 

We consider a cigar-shaped BEC tightly confined in the radial (x,y) plane, irradiated by a plane- 
wave laser as in [3] (see Fig. The laser polarization is along the long z-axis of the condensate to 
suppress collective ( "superradiant" ) Rayleigh scattering [201 or coherent atomic recoil lasing [SJ that 
are forbidden along the direction of polarization. The tight confinement along the propagation direction, 
together with the far off-resonance condition, enables us to treat the electromagnetic field within the 
Born approximation. The dipole-dipole potential induced by far off-resonance electromagnetic radiation 
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of intensity /, wave- vector Icl = fay (along the y-axis), and polarization e = z (along the z-axis) is [22| 

tfddM = J " ^ 4 [(! - 3 cos 2 (6>)) ( cos (far) + far sin (far) ) - sin 2 (e»)fc^r 2 cos (far) 1 cos(fay). (2) 

Here r is the interatomic axis, and 9 is the angle between r and the z-axis. The far-zone (far 
1) behavior of (J2J) along the z-axis is proportional to — sin(far)/(far) 2 . In terms of the condensate 
density n(r) = |i/>(r)| 2 , the mean-field energy functional accounting for atom-atom interactions is taken 
to be the sum H dd + H s [HQS 01, where Hdd = ^/ 2 ) / n ^ Udd ( r ~ n(r') d 3 r d 3 r' , and H s = 
(l/2)(4:TTah 2 /m) J n 2 (r)d 3 r is due to short-range (r -6 ) van der Waals interactions, which are described, 
as is usual, by a delta function pseudo-potential U s (r) — (Anah 2 /m)S(r). In momentum space Hdd takes 
the form, H dd = (1/2)(2tt)- 3 / U dd (k) n(k) n(-k) d 3 k, where U dd {\i) = / exp[-ik ■ r] U dd (r) d 3 r is the 
Fourier transform of the dipole-dipole potential J2Jl, the explicit form of which, for the laser propagation 
and polarization as in Fig.^ we have given in Eq. (4) of 

When the radial trapping is sufficiently tight so that the BEC is in its radial ground state we may 
adopt a cylindrically symmetric gaussian ansatz for the density profile whose width, w r , is 



3D (r) = Ninw 2 )-^) exp [-{x 2 + y 2 )/w 2 ] , (3) 



n 

where N is the total number of atoms and n(z), the ID axial density, is normalized to 1, but is un- 
specified for the time being. Denoting by n 3D (k) the Fourier transform of n 3D (r), we have n 3D (k) = 
Nn(k z )exp [— w 2 ,(k 2 + fc 2 )/4], in which n(fa) is the Fourier transform of the axial density n(z). This 
ansatz allows the principal value of the radial integration in H dd to be evaluated analytically so that the 
dipole-dipole energy reduces to a one dimensional functional along the axial (z) direction 

H dd = (N 2 /2) J n{z)n(z')U z dd (z - z') dz dz 1 = {N 2 /Ak) J n(k z )n(-k z )u[ d {k z ) dk z . (4) 

Defining the variables 

( = k 2 w 2 r /2 , £ = (k 2 - kt)w 2 r /2, (5) 
the one-dimensional (ID) axial potential that appears in H dd is 

u[ d (k z ) = I £^Q{t(w ri k z ),C(w r )}, (6) 

where Q[£(w r , fa), C(«v)] = ^ [-§ + F(£, 0] , and 

F(£, C) = 2^ exp(£ - C) £ ~ K{^i+i(0>- (7) 
j=o h 

%l{Ej[z]} is the real part of the generalized exponential integral [231 • However, the sum in the expression 
is rather unwieldy to use in calculations, so instead we will substitute it by the following asymptotic 
representation for the function F: 



M 

c+e' 
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This representation is very good for £ > and gives the correct asymptotics for £ — * (k z — > fcj,) and 
£ — > — C — ► 0). For — C < £ < the agreement is reasonably good (see Figure|5J|. 

Using the above expression, the Fourier Transform of the total (s-wave plus dipole-dipole) ID reduced 
interatomic potential is 

U t z ot (k z )=AE R a((k L w r )- 2 +lQ(w r ,k z )) (9) 

where En = h 2 k^/2m is the photon recoil energy of an atom and X is the dimensionless laser 'intensity' 
parameter 

It is emphasized that the radial degree of freedom is contained in J3J via the radius w r . Note that the 
Fourier Transform of the total potential is an even function: U^ ot (k z ) = U* ot (—k z ). 



3 Hartree-Fock-Bogoliubov-Popov method 

The Hartree-Fock-Bogoliubov (HFB) method is concisely described in ^2]. A reasonable approach within 
the HFB method is the Popov approximation JT] , which neglects the contribution of the above-condensate 
fraction to the anomalous correlation function. The advantage of the Popov approximation is that it gives 
a gapless spectrum of elementary excitations of a degenerate bosonic system in a wide temperature range 
(up to the critical temperature). In the present paper we apply the Popov approximation to a quasi- ID 
BEC with the laser-induced dipole-dipole interactions at finite temperature. Of course, for an infinitely 
long ID gas there is no condensate, even if this gas is non-ideal. This is because the depletion diverges 
at infinitely long wavelength. However, taking into account the finite length £ of the sample, we set an 
effective cutoff at a wavelength ~ £ and thus remove the divergence. Hence, a macroscopic population of 
the ground state becomes possible (in the case of an interacting gas). 

The ID ansatz for the atomic field operator is &(x,y,z,t) = t/)(z, t) exp[— (x 2 + y 2 )/(2w 2 )]/(y/Trw r ). 
Then the ID Hamiltonian reads as 

H = J dz^(z) f-L-^j^ + l JdzJ dz'$Hz')tHz)Ul t (z-z')$(z)iP(z'), (11) 

where U^ ot (z) is the inverse Fourier Transform of U^ ot {k z ) defined by Eq. (JHJl. Decomposing the atomic 
field operator in plane waves ip(z) = ^2 k £~ x l 2 exp(ikz)a k we rewrite Eq. (Ill|) in the momentum space: 

7j2^,2 y - 

= X, l^f 3 * 6 * + 2£ ^ U toMal +q a\,_ q a k ,a u . (12) 

k k,k' ,q 

The Hamiltonian of Eq. i|12|) applies to a system with a fixed number of particles. For simplicity we will 
work within the grand canonical ensemble and use the chemical potential, fi, to fix the mean particle 
number N . Then the Heisenberg equation of motion for the atomic annihilation operators d k , is given 
by the commutator of d k with H — nJ2 k ®kdk, an d reads 



d ( ?i 2 k 2 \ 1 
ih—d k = ( - n) a k + j^U? ot (q)a ] k , +q a k >a k+q . (13) 

k' ,q 
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Since the state with zero momentum is macroscopically populated, we invoke broken symmetry arguments 
[5] and replace clq by a c-number y/N c exp(itp), where N c is the number of atoms in the condensate. The 
phase ip is not an observable quantity, and its choice is arbitrary, thus we set <p = 0. The linear (ID) 
condensate density is defined as n c = N c /£. Atoms with k ^ comprise the above-condensate fraction 
(i.e., depletion) with the linear density 



i 

k^O 

which, together with n c , yields the total linear density 

n = n c + n a = N/£. 

Finally, in the Popov approximation the chemical potential is given by the expression 

fi = U^ t (0)(n c + n a )+W(0), 

where 

W(k) = - ^t(fc') 



(14) 



k'^-k 



a\, +k a k '+k 



(15) 
(16) 



For the modes with fc^Owe obtain, from Eq. (|13|l 



in>—ak = 
at 



2l2 



2m 



n c a k 



U^(0)n a + W(k)] a k + Uf ot (k)n c al k . (17) 



a k ak ) S k > k . Also note that the anomalous 



When deriving Eq. I|17|) . we take into account that yi\,a k i 

correlation functions for the above-condensate operators, such as (a k a k '}, are neglected in Eqs. i(ToHT7|) . 
as required by the Popov approximation. 

We now implement the standard Bogoliubov transformation |S] 

a k = u k exp(-iuj k t)b k - v k exp(iw fe t)fe^ fc , 

where the new (quasiparticle) operators b\, b k obey the usual bosonic commutation relations. The solution 
of the generalized eigenvalue problem gives the dispersion relation 



where 



huj k = JT(k) T(k)+2U t % t (k)n c , 



T(k) = %^ + W(k)-W(0). 



2m 



The transformation coefficients are 



u k 



; r(/,-i + P i ; ll (/,')r^ + l 



2hu) k 



v k 



< T{k) + U* ot {k)n a 1 
2hu; k 2 ' 



(18) 



(19) 



(20) 
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Note that if U* ot (k) were fc-independent (as it is in the case of absence of the external laser radiation, 
1 = 0) then T(k) would coincide with the free atom kinetic energy, and the well-known dispersion 
relation for a BEC with short-range interactions in the Popov approximation |121 ITT] would hold. At a 
finite temperature O (measured in energy units) we have |12| 



The Hartree-Fock-Bogoliubov-Popov equations should be solved self-consistently to obtain the con- 
densate and above-condensate occupation numbers, under the constraint that the total density is fixed, 
ri = n c + n a . This involves an iterative procedure |12| : 

1) As a starting point we choose the occupation numbers given by the {uk,Vk} coefficients found from 
the basic Bogoliubov scheme [5] (which does not treat the depletion self-consistently). 

2) The occupation numbers are used to calculate W(k) [Eq. I|16l) ] and hence the dispersion relation hu>k 
[Eq. (PJ)]. 

3) The occupation numbers are summed to give the new total density n', which is used to renormalize the 
occupation numbers of the condensate and above condensate modes, in order to keep the total density 
fixed. 

4) The dispersion relation, incorporating the new occupation numbers, is then used to recalculate the 
{uk,Vk} coefficients [Eq. J2UJ]- 

5) Repeat above steps 2-4 until convergence is achieved. 

4 Results and Discussion 

We have solved the HFB equations with the Popov approximation, as described above for the basic setup 
shown in Figure^ In Figures |3HH we present the results for two sets of parameters: 
Case 1: k L w r = 2.4, na = 0.66, £/w r = 74, Q/E R = 0.68, J = 1.12; 
Case 2: k L w r = 8.0, na = 4.7, £/w r = 28, Q/E K = 2.3, J = 0.28. 
For example, for 87 Rb these sets of parameters mean explicitly: 

Case 1 : The condensate radius w r is set at 0.3 /jm, the 3D density at the trap centre is equal to 
4 • 10 14 cm~ 3 , I = 22 um, 9 = 120 nK. The laser intensity can be found from the relation I/A 2 = 
0.26 W/(cm • GHz) , where A is the laser detuning from resonance. 

Case 2: w r = 1.0 /im, the peak 3D density is the same as in case 1, £ = 28 /im, 9 = 410 nK, 



Case 1 is representative of a BEC at the crossover between the ideal-gas and Thomas-Fermi limits with 
respect to the radial motion, while in case 2 the radial motion satisfies the conditions of the Thomas- Fermi 
limit. The parameters in cases 1 and 2 are chosen so that a significant 'roton' dip appears only in the 
more sophisticated Hartree-Fock-Bogoliubov-Popov method, but not in the simple Bogoliubov approach. 

In each of Figures EHU we have plotted three curves. The solid line gives the results of the Hartee- 
Fock-Bogoliubov-Popov calculation described in section [3] At convergence we found that in case 1 55 % 
of atoms remain in the condensate and in case 2 the condensate fraction is 82 %. The two other curves 
in the figures give the results obtained within the simple Bogoliubov approach and differ in the density 
that appears in the Bogoliubov dispersion relation. The variant (a) uses as the density in the Bogoliubov 
dispersion relation the condensate density n c obtained in the numerical calculations. The variant (b) uses 
in the same context as the density the total density of the system, n = n c + n a . Similarly, we calculate 




(21) 



I/A 2 = 0.064 W/(cm • GHz) 2 . 
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the depletion within the simple Bogoliubov approach for the variants (a) and (b). The position of the 
dip centre remains unchanged with temperature. 

Figures |21 and [SI show that when the depletion is taken into account self-consistently, the 'roton' dip 
in the excitation spectrum can actually become deeper. This means that the 'roton' dip persists at 
finite temperatures (of course, below the condensation temperature). Thus the sample heating due to 
spontaneous scattering of laser photons need not prevent observation of the 'roton' dip in the spectrum 
of the elementary excitations. Contrary to what one might expect, finite temperature effects not only 
do not wash out or diminish these intermode correlations, but rather serve to enhance them. This is 
surprising, considering the noisy character of finite-temperature density fluctuations. This might in part 
be because the dipole-dipole potential has a momentum dependence such that it is more effective at finite 
interatomic momenta (see Fig. 

One may ask: what happens if the laser intensity is increased so that the roton dip passes through the 
zero-energy axis, so that we obtain imaginary excitation frequencies? Two scenarios are then possible. 
The first one is that this would create a dynamical instability similar to that found in a BEC of atoms 
with negative scattering length |24j . The instability would grow exponentially, finally leaving us with 
a totally depleted BEC. Another scenario is the following: when the laser intensity reaches the value 
at which the roton minimum touches the zero-energy axis we may access a new ground state that is 
periodically modulated along the z-axis: a supersolid. An analogous behaviour had been predicted for 
superfluid liquid helium passing through a pipe [251. The formation of a supersolid state in a BEC 
with laser-induced dipole-dipole interactions has recently been studied j^HI by numerically integrating 
the time-dependent generalized Gross-Pitaevskii equation describing a cigar-shaped BEC irradiated by 
a circularly polarized laser (in contrast to the linear case considered in the present paper), whose wave 
vector is along the major axis, z, of the cigar. This new supersolid state has to be stable against 
small perturbations. One can choose between these two scenarios by numerical investigation of the BEC 
dynamics in the case of gradually-increasing laser intensity, based on a time-dependent generalization of 
the Popov approximation. Such an investigation is beyond the scope of the present paper. 

Finally, we would like to briefly mention experimentally detectable signatures of the roton minimum. 
Whilst the most convincing method would certainly be to map out the dispersion formula using Bragg 
spectroscopy, as in the measurement by Steinhauer et al [5], Bragg spectroscopy is technically difficult 
and time consuming. A simpler method would be to perform a time-of-flight measurement after releasing 
the BEC from the trap in order to detect the enhanced number of atoms having momenta around the 
roton minimum. 

To summarize, we have presented a method allowing self-consistent calculation of the effects of laser- 
induced dipole-dipole interactions upon a BEC in the presence of significant condensate depletion. This 
is especially relevant to the experimental situation where the temperature cannot be significantly lower 
than the chemical potential. This may help us cope with the adverse effects of the Rayleigh scattering of 
laser light [2j, which tend to heat up the sample. 
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POLARIZATION 




Figure 1: The laser beam and condensate geometry. 




Figure 2: Comparison of the exact (dashed line) and approximate (solid line) expressions for F(£, £) for 
(1) £ = 2.9, the solid and dashed lines arc practically indistinguishable, (2) £ = 32, the difference between 
the two lines is small but visible. 
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Figure 3: The spectrum of the elementary excitations of a BEC with laser- induced dipole-dipole interac- 
tions, for the set of parameters of case (1), see text. The solid line corresponds to the exact numerical 
solution, the dashed and dot-dashed lines correspond to the Bogoliubov spectrum for the variants (a) 
and (b), respectively. The frequency unit is the recoil frequency = En/h. 




Figure 4: The depletion (above condensate fraction) of a BEC with laser-induced dipole-dipole interac- 
tions, for the set of parameters of case (1), see text. The solid line corresponds to the exact numerical 
solution, the dashed and dot-dashed lines correspond to the depletion calculated for the variants (a) and 
(b), respectively. Note the correspondence between the position of the peak in the mean occupation 
number and that of the 'roton' dip in Fig. (see Ref. 3 ). 
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Figure 6: The same as in Fig. 0]for the set of parameters of case (2), see text. 
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